Adapt3D.h

Go to the documentation of this file.

/*
 * Adapt3D.h
 *
 *  Created on: Oct 29, 2010
 *      Author: David Berrios
 */

#ifndef ADAPT3D_H_
#define ADAPT3D_H_
#include <string>
#include "Cell3D.h"
#include "Model.h"
#include "Constants.h"
#include "Point3f.h"
#include <vector>
#include <boost/unordered_map.hpp>
#include <boost/multi_array.hpp>
#include <string>

#define NNODE_ADAPT3D 4
#define NVARS_ADAPT3D 12
#define NDIMN_ADAPT3D 3
#define INNER_RADIUS 1.0
#define OUTER_RADIUS 30.0

#define nx_sg 100
#define ny_sg 100
#define nz_sg 100
#define nx_b 2
#define ny_b 2
#define nz_b 2

//#define CARTESIAN_S_GRID
#define SPHERICAL_S_GRID
#define DELAUNEY_SEARCH
#define DELAUNEY_ITER_MAX 1000                                          /* new 11-16-11 */


namespace ccmc
{

    struct SmartGridSearchValues
    {
        float xl_sg,xr_sg,yl_sg,yr_sg,zl_sg,zr_sg;
        float xl_gr,xr_gr,yl_gr,yr_gr,zl_gr,zr_gr;
        float dx_sg,dy_sg,dz_sg;
        std::vector<int> * indx;
        std::vector<int> * esup1;
        std::vector<int> * esup2;
        boost::multi_array<int, 3> nelems_in_cell;//(boost::extents[nz_sg][ny_sg][nx_sg]);
        boost::multi_array<int, 3> nnodes_in_cell;//(boost::extents[nz_sg][ny_sg][nx_sg]);
        boost::multi_array<int, 3> start_index;//(boost::extents[nz_sg][ny_sg][nx_sg]);
        boost::multi_array<int, 3> start_index_nodes;//(boost::extents[nz_sg][ny_sg][nx_sg]);
        boost::multi_array<int, 3> end_index;//(boost::extents[nz_sg][ny_sg][nx_sg]);
        boost::multi_array<int, 3> end_index_nodes;//(boost::extents[nz_sg][ny_sg][nx_sg]);
        std::vector<int> * indx_nodes;
        //int still_in_same_element;
        std::vector<int> * delauney_search_iteration_profile;
        std::vector<int> * facing_elements;
        //int outside_grid;
        int last_element_found;
        //BoundingBox * parent;
    };
    class Adapt3D: public Model
    {

        public:
            Adapt3D();
            long open(const std::string& filename);

            Interpolator * createNewInterpolator();
            SmartGridSearchValues * getSmartGridSearchValues();
            //float * getModifiedCoords();
            const std::vector<std::string> getLoadedVariables();

            virtual ~Adapt3D();

            /*template<typename T>
            friend int binary_search(const std::vector<T>& vec, unsigned int start, unsigned int end, const T& key);*/

        protected:
            void initializeConversionFactorsToSI();
            void initializeSIUnits();

        private:
            int ndimn;
            int nnode;
            int nvars;
            int grid_reg_no, npoin, nelem, nboun, nconi;
            long GRID_REG_NO_cdf_num, NPOIN_cdf_num, NELEM_cdf_num, NDIMN_cdf_num;
            long NBOUN_cdf_num, NCONI_cdf_num;

            void smartSearchSetup();
            bool setupSearchUnstructuredGrid();
            void setupOctreeGrid();
            std::vector<Point3f> vertices;
            SmartGridSearchValues smartSearchValues;
            std::vector<float> *coord_modified;
            std::vector<int> * intmat_modified;
            int index_2d_to_1d( int i1, int i2, int n);

            void locate_facing_elements();






    };
}

#endif /* ADAPT3D_H_ */